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An analysis method of lifetime, polarization and spectrally fil- 
tered fluorescence correlation spectroscopy, referred to as fil- 
tered FCS (fFCS), is introduced. It uses, but is not limited to, 
multiparameter fluorescence detection to differentiate be- 
tween molecular species with respect to their fluorescence life- 
time, polarization and spectral information. Like the recently 
introduced fluorescence lifetime correlation spectroscopy 
(FLCS) [Chem. Phys. Lett. 2002, 353, 439-445], fFCS is based on 
pulsed laser excitation. However, it uses the species-specific 
polarization and spectrally resolved fluorescence decays to 
generate filters. We determined the most efficient method to 
generate global filters taking into account the anisotropy infor- 
mation. Thus, fFCS is able to distinguish species, even if they 
have very close or the same fluorescence lifetime, given differ- 
ences in other fluorescence parameters. fFCS can be applied as 



a tool to compute species-specific auto- (SACF) and cross- cor- 
relation (SCCF) functions from a mixture of different species for 
accurate and quantitative analysis of their concentration, diffu- 
sion and kinetic properties. The computed correlation curves 
are also free from artifacts caused by unspecific background 
signal. We tested this methodology by simulating the extreme 
case of ligand-receptor binding processes monitored only by 
differences in fluorescence anisotropy. Furthermore, we apply 
fFCS to an experimental single-molecule FRET study of an 
open-to-closed conformational transition of the protein Syntax- 
in-1. In conclusion, fFCS and the global analysis of the SACFs 
and SCCF is a key tool to investigate binding processes and 
conformational dynamics of biomolecules in a nanosecond-to- 
millisecond time range as well as to unravel the involved mo- 
lecular states. 



1. Introduction 

Fluorescence correlation spectroscopy (FCS) is frequently used 
to monitor molecular interactions (binding and unbinding), 
chemical reaction kinetics and diffusion constants of fluores- 
cent molecules. 11 " 91 FCS, developed in the early 1970s, 11,101 im- 
proved significantly in the 1990s due to better optics and exci- 
tation sources. 111,121 Regardless of its ultrasensitive detection, 
complexities arise when taking into account a mixture of fluo- 
rescent molecular species with different brightness and differ- 
ent diffusion constants, in particular when the diffusion con- 
stants are similar. 1131 Furthermore, scatter background is known 
to reduce the accuracy of determination of the true number of 
molecules in the observation volume. 1141 Several ideas to over- 
come these problems have been implemented in literature. 
One approach is to combine pulsed excitation and time-re- 
solved fluorescence detection with FCS and time-gate the fluo- 
rescence signal in the nanosecond range. This approach can 
be used to suppress light scattering, but also it can help to 
separate a dye with a short fluorescence lifetime with respect 
to a dye with a long fluorescence lifetime. 115 " 211 However, one 
major disadvantage of the time-gate approach is that all pho- 
tons arriving outside the time-gate are discarded in the analy- 
sis and becomes too subjective in definition. Another disad- 
vantage is that it cannot fully remove unwanted signal be- 
cause of "crosstalk" into time gate. 

Recently, a different strategy using different fluorescence 
lifetimes of dyes was proposed and named fluorescence life- 
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time correlation spectroscopy (FLCS). 122 " 251 In this case, two in- 
dependent FCS curves for two different dyes, with different 
lifetime decay, were extracted from a mixed solution without 
ignoring any detected photons. The authors describe the 
method as being "similar to recording a complete spectrum 
and deconvoluting it into the different spectral contributions 
from the various emitting species". One of the advantages of 
this method is that there is no lost signal. In addition, one can 
discriminate detector after pulsing from the fluorescence 
signal. 1231 But one of the most important advantages is that 
one can obtain FCS curves for different dyes. This method can 
be generalized for more than two dyes, the only requirement 
is that they have sufficiently distinct fluorescence lifetimes. 1221 
Ruttinger et al. 1261 demonstrated that diffusion times differing 
only 25% from each other can be resolved by FLCS. Various 
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applications of FLCS and fluorescence lifetime crosscorrelation 
spectroscopy (FLCCS) demonstrated the usefulness of this 
technique for characterizing species-specific diffusion times 



and biomolecular binding processes in vitro, 1 



in live cells 1291 



and at a light absorbing interface,™ which employed fluores- 
cence lifetime tuning to study diffusion of BODIPY-tail-labeled 
lipid in distinct lipid bilayer environments. Moreover, FLCS was 
applied to study the protonation kinetics of fluorescent dyes in 
aqueous buffers. 1311 

The requirement of FLCS, that fluorescence decays are suffi- 
ciently distinct due to different fluorescence lifetimes, becomes 
one of its major limitations. In practice — from our experi- 
ence — "sufficiently different decay times" means at least 1.0- 
1.5 ns difference for efficient separation. This corresponds to 
20-30% change in lifetimes of fluorescent dyes, commonly 
ranging from 1 to 5 ns. Therefore, this requirement significantly 
restricts dye selection, the number of species, which can be 
studied in parallel, and the type of problems that could be 
studied. 

The methodology proposed in this work, fFCS, benefits from 
multiparameter fluorescence detection (MFD). 116,32,331 Using 
single-color excitation and MFD, the collected signal is split 
into several detection channels at different spectral windows 
considering their parallel and perpendicular polarizations rela- 
tive to linearly polarized excitation laser beam. In addition to 
fluorescence lifetime and time resolved anisotropy we can also 
utilize the spectral differences of species as additional parame- 
ter to further improve the distinction between molecular spe- 
cies in a mixture. Two or three dimensional filtering (lifetime 
and anisotropy, lifetime and color, or lifetime, anisotropy and 
color) makes this method more precise and selective, for exam- 
ple, 16 species have been studied by Widengren at al. [331 In par- 
ticular, it becomes possible to "deconvolute" FCS curves for 
the contribution of different molecular species even when they 
have very similar or equal lifetimes. The only requirement is 
that they differ sufficiently in one parameter or in a combina- 
tion of parameters if their individual differences are small. Here 
we demonstrate the case in which a difference of 0.2 ns in ro- 
tational correlation times is enough to separate species even 
when their fluorescence lifetimes are identical. 

We tested fFCS for two cases: 1) To characterize the accuracy 
of fFCS, we simulated the dynamic binding of two biomole- 
cules (e.g. labeled ligand and unlabeled receptor). From fFCS, 
we extracted what we call the species auto-correlation func- 
tion (SACF) and species cross-correlation function (SCCF) to 
quantitatively resolve species fractions in dynamic equilibrium 
and recover rate constants with relative error less than 2% in 
average. The separation of species is possible even in the ex- 
treme case where the two biomolecules have the same fluo- 
rescence lifetime, but differ by rotational correlation time. Time 
resolution of sub-microseconds can be achieved. 2) Further- 
more, we applied this methodology to recent experimental 
data 1341 to extract kinetic information about conformational 
transition between the open and closed state of the protein 
Syntaxinl with high temporal resolution. In this case, separa- 
tion of conformers is maximized by the differences in lifetime 
and anisotropy using single molecule Forster Resonance 



Energy Transfer (smFRET) data. Finally, the requirements for an 
optimal experimental realization of fFCS are discussed. 

Combining FRET and FCS has been attempted before; 135 " 371 
however, there are multiple artifacts that complicate the analy- 
sis and limit its use. For example, it is almost impossible to 
avoid donor-only (D-only) labeled species in mixture with 
donor and acceptor (DA) labeled species. On top of this, one 
can have photophysical processes, such as triplet formation, 
cis-trans isomerization, and donor quenching. Therefore, the 
detection of small FRET efficiency changes due to conforma- 
tional transitions in FRET-FCS is a very difficult task due to 
poor contrast. Nevertheless, with fFCS in principle these com- 
plexities can be solved. Additionally, one could apply fFCS in 
other variations of single molecule FRET experiments, such as 
pulsed excitation with several lasers referred to as nsALEX 1381 
and PIE 1391 (or see article by Kudryavtsev et al. published in this 
same issue of ChemPhysChem m ), and to assure that the ac- 
ceptor dye is functional and in this way increase selectivity and 
confidence that the monitored species are relevant to the 
study. 

Another advantage of fFCS is to accurately separate pure 
fluorescence from background signal, in particular when meas- 
uring pM concentrations in single-molecule experiments. We 
show in this report, that fFCS is successfully applied to obtain 
correlation curves with correct amplitudes from single mole- 
cule data. 

In all, fFCS allows one to extract all possible auto- and cross- 
correlation functions that could be used to study a wide varie- 
ty of processes including dynamic binding and conformational 
dynamics of proteins. 



2. Results 

2.1. Finding Optimal Filters 

The auto-correlation function allows direct assessment of the 
diffusion constant D. If only diffusion of a single species is con- 
sidered, the amplitude at zero time of the autocorrelation 
function G(0), allows one to determine the mean number of 
molecules N in the detection volume, V det , or the concentra- 
tion, c, if the parameters of detection volume are known 
[Eq. (1)]: 
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Unfortunately, the situation becomes significantly more com- 
plicated if more than one molecular species are simultaneously 
present in solution. Considering a mixture of n species, with 
corresponding brightnesses Q w , diffusion constants D (/) and 
fractions x M , /'= 7,..., n, it is most convenient to define the au- 
tocorrelation function as 1401 [Eq. (2)]: 
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where c H concentration, Q (,) brightness and x (,) 



fraction 



of species /', and N corresponds now to the total number of 
molecules. The normalized diffusion term of species / [i.e. 
G di ff (u) = 1] is given by Equation (3): 



ll p (/) + V ) = ^| lp (0 + £ ip (0 = £ p W = 1 



(4) 



where the parallel channel probability is [Eq. (5)]: 



(3) llpM 



This model assumes a three-dimensional (3D) Gaussian-shaped 
volume element with spatial distribution of the detection 
probabilities w(x,y,z) = exp(— 2(x 2 + y 2 )/a>o) exp(— 2z 2 /z 2 ). 
The 1/e 2 radii in x, y and in z directions are denoted by a> 0 and 
z 0 , respectively. For one-photon excitation, the characteristic 
diffusion time tf can be used to estimate the diffusion coeffi- 
cient 0<° by tf = « 2 /4D« 

A simple relation like Equation (1) cannot be found for multi- 
ple species. Even if one knows the brightness Q w for each spe- 
cies, the diffusion coefficients Distill have to be significantly 
distinct for successful extraction of c (/) values, highlighting the 
need of a proper methodology that could differentiate species 
in a mixture. In ref. [22] it was proposed to use the differences 
in fluorescence lifetime to separate molecular species by time- 
correlated single-photon counting (TCSPC). 141,421 For simplicity, 
we will describe how to generate optimal filters for separating 
two molecular species with mono color detection split for po- 
larization. However, fFCS can be easily extended to account for 
spectral differences. In addition, fFCS can in principle be used 
to separate more than two species. 

To fully use the anisotropy information a relationship be- 
tween the total intensities in the parallel and perpendicular de- 
tection channels is kept. Therefore the methodology to find 
the proper filters by minimizing the relative errors proposed in 
ref. [22] needs to be optimized. Having this in mind, there are 
in principle three possible scenarios to generate filters: Scenar- 
io 1: single detector, single filter {' d ~ s f). Scenario 2: two detec- 
tors, independent filters ( 2d ~'"f). Scenario 3: two detectors, 
global filter ( 2d ~ 9 'f). We found that scenario 3 is most efficient in 
separating the species. It is described in more detail below. 
The theory behind scenario 1 was described in ref. [22], and in 
Section 4.1 we present the case for scenario 2. 

In the case of MFD, where the situation is not as trivial as in 
the case of one-channel detection, the fluorescence signal 
from the mixture is divided into its parallel and perpendicular 
components and into two spectral ranges. Each registered 
photon emitted by fth species is detected with probabilities 



p w in either parallel or perpendicular detection channel, 

L pf exist due to specif- 



Y\ , 

and certain relations between "pj' 1 and 
ic anisotropy decays. If these relationships are disregarded (as 
done in scenarios 1 and 2), the anisotropy information is lost. 
To solve this problem, in scenario 3 we stack the TCSPC histo- 
grams of the perpendicular over the parallel detection chan- 
nels and apply a global normalization to the stacked TCSPC 
histogram. As result, the corresponding probabilities of each 
fth species are obtained as [Eq. (4)]: 
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and similarly, the perpendicular channel probability is [Eq. (6)]: 

(6) 



Thus, "pj' 1 and x pf represent the conditional probabilities to 
register a photon in the ;'th bin of parallel or perpendicular 
TCSPC channels, provided that photon is emitted by fth spe- 
cies. To differentiate from what is done in Section 4.1, here we 
define pf as the stacked conditional probability of length 2L 
This is done to maintain the anisotropy information. To ac- 
count for spectral differences (color detection) the histograms 
are stacked accordingly. 

The total number of registered photons W is the sum of 
photon numbers emitted by all species in mixture solution 
[Eq. (7)]: 
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where w (t> is the photon counts of species /. Instead of two 
measured decay histograms ("/-/, and x Hy) of the two species in 
the mixture we consider one bimodal decay histogram with 
total number of the stacked TCSPC channels 21 and the condi- 
tional probability distribution pf for the fth species. This bimo- 
dal decay distribution can be defined as [Eq. (8)]: 
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Where [Eq. (9)]: 
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Using the global definition in Equation (8), two bimodal filters 
rj , /=1, 2, will be generated with the property [Eq. (10)]: 



(10) 



where the brackets denote averaging over long time of meas- 
urements. Filters f-' ] are obtained by minimizing relative 
errors 1431 for parallel and perpendicular detection channels si- 
multaneously [Eq. (11)]: 
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The idea of one bimodal decay histogram Hj and two condi- 
tional probability distributions pf, one for each species with 
total number of stacked TCSPC channels 2L, makes it possible 
to calculate what we refer to as the species auto- (SACF) and 
cross-correlation (SCCF) function in a similar fashion to 
ref. [22]. However, we need to transform our raw data streams 
into a modified format where the TCSPC channels of both de- 
tectors are stacked in a single array with length 2L In this way, 
the SCCF G (, ' m) (f c ) between species i and m is [Eq. (12)]: 



G (/J7 "(r/ 



(12) 



where S ; (f) is the signal in the y'th stacked TCSPC channel array 
of the total signal at measurement time f and Sj(t + t c ) is the 
signal in the stacked TCSPC channel y'th at measurement time 
f + f c . The SACF is defined for the case of i=m while SCCF cor- 
responds to different species i^m. In both cases, the ortho- 
normality relationship [Eq. (13)]: 



7=1 



/ = m, 
i^=m. 



(13) 



is satisfied. Then, the species auto- and cross-correlation func- 
tions G (,m) (t c ) averaged over an infinite number of measure- 
ments or over sufficiently long measurement time (>t c ) are 
[Eq.(14)]: 



G<'»(f c ) 



(F")(t) •F (m) (r + Q) 
(F«(t))-(F^(t + t c )) 



(14) 



where F (,) (f) and F (m) (f + f c ) are pure fluorescence signals from 
the molecular species of rth and mth type. In a two-state 
system, G (,m) (f c ) and G (mi) (f c ) are equal. Thus, since the aniso- 
tropy differences are counted in addition to lifetimes, it be- 
comes possible to highlight any dynamic process between two 
species in mixture solution if they differ in rotational correla- 
tion time or/and lifetime. In contrast to standard correlation 



curves the amplitude of dynamic term per molecule in species 
cross correlation function is equal to —1, like in an antibunch- 
ing term that will be discussed in Section 2.2.1. One can under- 
stand the shape of SCCF as a probability distribution of inter- 
conversion between species: for very short f c this probability is 
nearly equal to zero (SCCF starts at baseline level) and gets 
higher for larger f c proportional to relaxation time of dynamics. 
Finally, everything is limited by diffusion time and probability 
to observe a transition is dropping to baseline level. If dynam- 
ics between species is missing then SCCF is showing no addi- 
tional correlation amplitude above baseline. 

In comparison to FCS, the species auto-correlation function 
(SACF) as defined in Equation (12) with i = m assumes that 
other species are essentially non-fluorescent. In particular, it 
can be considered a background/scatter free autocorrelation 
function if a scatter filter was included. Then, the amplitude of 
the SACF is inversely related to the true number of molecules 
with the specific lifetime and polarization characteristics given 
by the filters. And in the same way as in standard FCS, at small 
lag times species are highly correlated and for longer correla- 
tion times everything is limited by diffusion showing a decrease 
of amplitude. 

Let us consider an extreme case where two molecular spe- 
cies are simulated for typical conditions in single-molecule ex- 
periment: both species have the same fluorescence lifetime 



(T, 



= 4.0 ns), same brightness Q (1) = Q' 



fD. 



differ only by rotational correlation times (p (1) = 0.1 ns " i2) 



150 kHz, but 

P° 

0.3 ns). Also, both species have the same diffusion time 
( f D = f D = 0-25 rns). The average number of simulated mole- 
cules in the observation volume was N = 0.399 (A/, =0.133 
and N 2 = 0.266). Details on the Brownian dynamics simulator 
can be found in Section 4.2. The chosen brightnesses are typi- 
cal for FCS using confocal microscopy. 144 " 461 By applying addi- 
tives like triplet or radical quenchers the signal can be further 
increased. 1471 

This plausible experimental situation is simulated to empha- 
size that only differences in anisotropy can be used to separate 
species. The same simulation data set is used to determine the 
SACFs of the three different scenarios of filter generation. In all 
cases, exactly the same number of photons is correlated. The 
detection is always considered in two detection channels (par- 
allel and perpendicular). Figure 1 shows the comparison of cal- 
culated SACF by Equation (12) curves for the three difference 
scenarios. In all panels solid blue and red lines correspond to 
the SACF for species 1 and 2 respectively. On top the modeled 
correlation curves are shown as dashed lines. 

Scenario 1: Single detector, single filter {' d ~ s fj. In this case, 
Figure 1 A, there is no split by polarization, although polarized 
excitation was used. To mimic a single detector from the same 
simulation, we assumed that the total signal is the sum of the 



two simulated channels h, 



"Hi 



Hi. The same 



hj decay pattern is used for both detection channels (parallel 
and perpendicular) and correspondingly the same filters are 
applied for parallel and perpendicular detection channels. The 
conditional probabilities "pj 1 ' = ± pj' ) are obtained via Equa- 
tions (4-6) using "fy and Filters are generated by simulta- 
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neously minimizing relative errors in both detection channels 
[Eq. (11)]. 

Scenario 2: Two detectors, independent filters ( 2d ~"Y). Fig- 
ure 1 B shows the SACFs using independent filter generation 
for each detection channel and stacking them afterwards. The 
conditional probabilities "pj 1 ' and ± pj' ) are obtained via Equa- 
tion (27). n fj and x ly'' filters are generated by simultaneously 
and independently minimizing relative errors in both detection 
channels [Eq. (29)]. 



Figure 1. Comparison of SACF obtained by using three possible filter-gener- 
ation scenarios. The species in the simulated mixture have the same lifetime 
but different rotational correlation times. Both species have the same bright- 
ness Q m = Q a) = 150 kHz, and identical diffusion times f d nl = f d l21 = 0.25 ms. 
The number of molecules per species are N'" = 0.133, N a> = 0.266, respective- 
ly. Dark counts = 0.2 kHz, and scatter (IRF) = 2.5 kHz rates were considered. 
The lifetime of both species is t g <1> = t g (2) =4.0 ns, and the respective rota- 
tional correlation times are p l1l = 0.l ns, p (21 = 0.3 ns. The total simulated ex- 
periment duration is 2479 s. The SACF for each species are shown in solid 
blue and red lines for species 1 and 2 respectively. For comparison the simu- 
lated correlation function of each species are shown as dashed lines in dark 
blue and wine for species 1 and 2. A) The calculated SACF by Equation (12) 
with no split in polarization (scenario 1: assuming "fy = ± /i J = "H, + ± H j as 
a case imitating single detection channel experiment). The same decay pat- 
tern is used for both detection channels and correspondingly the same fil- 
ters are applied for parallel and perpendicular detection channels to mimic 
the one detector case. The correlations show poor statistics with some 
degree of separation between species. B) scenario 2: SACF from Equa- 
tion (12) using two detectors, independent filter generation for each detec- 
tion channel and stacking them afterwards. The separation is better com- 
pared to the single detector case. C) SACF of scenario 3 using Equation (12) 
for two detectors and global filter generation. The separation of species is 
obvious to the eye, and the estimated error in the recovered parameters is 
within 2%. A detail error analysis of recovered parameters is done in Sec- 
tion 2.2. 

Scenario 3: Two detectors, global filter ( 2d ~ gl f). When using 
global filter generation the conditional probabilities "pj' 1 and 
x pj''are obtained via Equations (4-6). "/p and L fj filters are 
generated by simultaneously minimizing relative errors in both 
detection channels [Eq. (11)]. 

The case imitating one detection channel (Figure 1 A) shows 
poor statistics with some degree of separation between spe- 
cies because of polarized excitation. The separation gets better 
compared to the single detector imitating case for SACF from 
two detectors using independent filter generation (Figure 1 B). 
Compared to Figures 1 A, B, the contrast in Figure 1 C is greater 
and a clear separation of species is possible, if the anisotropy 
is properly considered. This is only possible when polarized de- 
tection, splitting the fluorescence signal into parallel and per- 
pendicular detection channels, is used like in single molecule 
MFD setups. In all following results scenario 3 of generating fil- 
ters is used. 

It is worth mentioning that in all cases, the background pho- 
tons were considered as a third species. The normalized TCSPC 
histograms {pf) for each species and their corresponding fil- 
ters according to the three different scenarios are shown in 
Figure 9. The use of the background filter is an additional stan- 
dard methodology used for correctly recovering the amplitude 
of the correlation and was used in all simulated and experi- 
mental data analyses. 

2.2. Results for Polarization-Resolved Simulations 

2.2.1. Static/Dynamic Binding Equilibria: Species Auto- and 
Cross-Correlation Function Separate Fractions 

To test the implementation of fFCS, we generated different si- 
mulated experiments using a Brownian dynamics approach 
(see Section 4.2). All fluorescence spectroscopic parameters of 
molecules freely diffusing through a confocal observation 
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volume were generated. First, we simulated a mixture of mo- 
lecular species which differ only by their rotational correlation 
times. This example could represent "static" mixture of free 
ligand molecules and static complexes (e.g. free labeled DNA 
and labeled DNA bound to an unlabeled protein as in the case 
of a binding experiment where the diffusion coefficient does 
not change significantly). Second, we simulated a mixture of 
free ligands and complexes, but now we consider a dynamic 
interconversion between these two species, which may repre- 
sent a kinetic exchange between a bound and free DNA. In 
the following we discuss the static case where the exchange 
relaxation time is much larger than the characteristic diffusion 
times (f R §>f d ) and the dynamic case with f R smaller or compa- 
rable to f d . A possible application of this example is to study 
conformational dynamics of biomolecules using a single dye, 
given that there is a significantly large change in its rotational 
correlation time. 

For the cases of static and dynamic binding equilibria Equa- 
tion (15), we have simulated mixtures of two species, free 
ligand, L f , and the complex, (LR), that is, ligand bound by the 
receptor R [Eq. (15)]: 



J2 

ra 
n 
o 



0.1 



0.01 



1E-3, 



1E-4 



IRF 

free ligand (p(Z_ f ) = 2.71 ns) 

bound DNA (p(LR) = 5.18 ns) 





Parallel channel 



Perpendicular channel 



256 512 768 

Stacked TCSPC channel # / 32.55 ps/ch 



1024 



(LR) 



The relative concentration of re- 
ceptor to ligand is considered 
high that the process can be 
modeled as a pseudo-first-order 
binding process with the rate 
constant of association k h2 = 
k a -[R] and rate constant of disso- 
ciation k v . The species fractions 
X in the dynamic equilibrium are 
defined by Equation (16): 



^1,2 



and 



(16) 



The inverse of the sum of the 
rate constants corresponds to 
the characteristic relaxation time 
of binding equilibrium f R observ- 
able in FCS [Eq. (17)]: 



(17) 



0. 



For the static case k l2 = l 
In other words, there is no ex- 
change between species. We 
have used the same 4 ns life- 
time for both species but differ- 
ent characteristic rotational cor- 



relation 



times, 



(free-ligand) _ 



2.71ns and p (complex) = 5.18 ns, 



(15) 



Figure 2. Conditional probabilities of polarization resolved fluorescence 
decays used in simulations for static and dynamic binding equilibria of two 
species in buffer: free ligand (p lfreel| 9 and > = 2.7l ns) and complex 

(p (comp """ = 5.18ns). 



are similar to experimental values from real DNA-protein bind- 
ing experiments. Both species had the same diffusion time, 
f d = f d = 1 ms, and the same brightness, 
q(M _ q(ir) _ t oo kcpm, (in order to make the situation very 
extreme). The concentrations with an average number N of 
molecules in the focus were chosen in the simulations such 
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Figure 3. Burstwise analysis of single-molecule events in simulations of an SMD experiment. 2D histograms of 
fluorescence lifetime (r G ) distributions on the x axis and scatter-corrected anisotropy r scG on the y axis for mixtures 
of two species [free-labeled ligand (p**** ll9sn * = 2.71 ns) and a complex (p (compl,9 ' l = 5.18 ns)] in buffer. The overlaid 
red and blue curves show the Perrin equation r scG = r 0 /(l + r G /p) with a fundamental anisotropy r 0 = 0.375 and 
the species-specific rotational correlation times p. A) 50 by 50 percent (X (LR '=0.5) mixture of static species. B) 80 
by 20 percent (X< LR> = 0.2) mixture of dynamic species (/c 12 = 2000 s~ 1 ; fc 2/ , =8000 s -1 or f„ = 0.1 ms), C) 50 by 50 per- 
cent (X"- R1 = 0.5) mixture of dynamic species (fc u = k 2 , = 2000 s" 1 or f„ = 0.25 ms). D) 20 by 80 percent (X"- R, = 0.8) 
mixture of dynamic species {k, 2 = 8000 s~'; k 2 , = 2000 s~ 1 or f R = 0.1 ms). The 1 D lifetime r a distributions are over- 
laid by corresponding ones from 50 by 50 percent static mixture (green line). The ID r scG distributions (gray) are 
overlaid by distributions from 100% free ligand (blue line), 50 by 50 percent mixture of static species (green line) 
and 100% complex (red line). 
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that they were optimal for the two measurement 
types: A/ = 0.1 for correlation analyses and A/=0.02 
for single-molecule experiments. The generated data 
files where then filtered with the time and polariza- 
tion decays shown in Figure 2. 

For direct visual control it is useful to check the be- 
havior of simulated model systems by MFD at single 
molecule level applying the standard analysis tech- 
niques described in refs. [32,48,49]. In Figure 3 the 
results of a burstwise analysis are displayed in two-di- 
mensional frequency histograms of scatter corrected 
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The corresponding one-dimensional parameter histo- 
grams are given as projections. As one can easily see 
from Figure 3, the lifetime distributions are not af- 
fected by dynamics. However, the differences in ani- 
sotropy histograms clearly show the expected mixing 
of states within single burst. As references the 1D r scG 
distributions (gray bars) are overlaid by distributions 
obtained from reference conditions: 1)pure free 
ligand (blue line); 2) pure complex (red line); and 
3) mixture with static equilibrium and X^^O.S 
(green line). For a static equilibrium (Figure 3 A), the 
anisotropy of the mixture is nicely described by the 
sum of the individual components. Considering dy- 
namic binding equilibria with an increasing amount 
of complex (Figure 3 B-D), the maximum of experi- 
mental distribution (gray bars) continuously moves 
from a lower r scG value (free DNA) to a higher one 
(100% complex), which indicates the increasing 
steady-state fraction of the complex. In Figure 3C, we 
see the presence of dynamics by the fact that the 
gray anisotropy distribution is narrower than the cor- 
responding static equilibrium (green line). However, 
these observations do not allow us to answer the 
question of how fast the binding dynamics is. In other words, 
we need to get a reliable technique to distinguish the dynamic 
binding from static equilibrium. Something that highlights the 
interconversion rate is needed, and fFCS provides such infor- 
mation. 

The standard FCCS curves (signal in parallel channel versus 
one in perpendicular) calculated from raw simulated TCSPC 
data for 50/50% static and 50/50% dynamic equilibria are 
overlaid in Figure 4A. Two features become immediately ap- 
parent — correlation amplitudes are lower than expected and 
there is no bunching term sensing the dynamic exchange. In 
the absence of background G(0) = 10 should be obtained since 
we know the exact number of particles 0.1 in excitation 
volume for simulated data. However, the contribution of un- 
corrected background from buffer results in lower amplitude 
of curves (G(0)=7) and wrong A/ FCS = 0.143. Moreover, regular 
FCS cannot detect the dynamics in this type of samples, and 
the decay is sufficiently described by a diffusion term (single 
gray bar at 1 ms in Figure 4A). Thus, the absence of any differ- 
ence may be easily misinterpreted as absolute "identity" of the 
two studied samples. This is a typical example where classical 
FCS alone very often overlooks properties of samples, which 



C) 


. Specie 


s Auto Correlation 




Wki E 


<-■ 


uiliOtium (X ^^O S) 

nd — □— ligand 
igand ■■ fH ligand 




fll "\ 


— fix 


jjmplex fil complex 




! V \ 








0.t\ms 

V t 













O) 

0.1 

3 


R Species Cross Correlation 

"J&f ^ x [ia) =a.5 

jFp \ 1 ms 

dt dynamic 
' -"- X m =02 

]r \ fitX lLI, ^0.2 

y \ -o- A" (LR> =0.5 

A \ fitX (LRl =0.5 

6 1r \ -o- X^0.% 

0.25 ms \ fitx ttR, =o.s 

'<': '-..',.^\ _ . ^**<w_ 


m 1 .f .. 





0.01 0.1 1 10 
Correlation Time / ms 



0.1 1 10 100 
Correlation Time / ms 



Figure 4. Simulation of an FCS experiment for distinct mixtures of species [free ligand 
(p <freeLi9andl = 2.71 ns) and complex (p ,compl,9<l = 5.18 ns)] in buffer with W = 0.1. The data 
generated in the simulation correspond to a measurement time in a real experiment of 
approximately 12500 s. The geometric shape of the excitation focus was defined as 3D 
Gaussian with co 0 = 0.2 |im and z 0 = 0.6 |im. The background signal consists of dark 
counts = 0.2 kHz and scatter= 1.8 kHz. The filtered FCS curves were computed by using 
the filters generated from polarization-resolved decays shown in Figure 2. Comparison of 
FCCS, SACFs and fit curves for the following cases: A) Overlay of FCCS curves calculated 
from raw simulated data for 50/50% static (gray) and 50/50% dynamic mixtures (green). 
B) Overlay of SACFs for 50/50% mixtures: static and dynamic equilibrium, respectively, 
with (k, 2 = k 2 , = 2000 s _1 or f R = 0.25 ms). C) Overlay of SACFs for 20/80% mixtures: static 
and dynamic equilibrium, respectively, with (fr 12 = 8000 s -1 , fr 21 = 2000 s~ 1 or f R = 0.1 ms). 
D) Overlay of SCCFs for two cases: i) 50/50% mixtures: static (gray curve) and dynamic 
(green curve, /(, 2 = /c 2 , = 2000 s~' or f R = 0.25 ms) equilibrium (fit results to f R = 0.246 ms); 
ii) 80/20% mixture in a dynamic equilibrium (pink curve, /( 12 = 2000 s -1 , k v =8000 s -1 or 
f B =0.1 ms). The fit gave f R = 0.096 ms; iii) 20/80% mixture in a dynamic equilibrium (blue 
curve, /( 12 = 8000 s~\ k 2 , =2000 s~' or f R = 0.1 ms). The fit results in f R = 0.099 ms. 



have been identified as heterogeneous by MFD (Figure 3). In 
Figure 4 B, we demonstrate that fFCS is free of these uncertain- 
ties. Applying the corresponding filter sets for three compo- 
nents in Figure 2 [background (IRF), species 1 (free ligand) and 
species 2 (complex)] obtained by global error minimization 
[Eq. (11)], the SACFs were computed using Equation (12). In 
Figures 4B,C, the SACFs for static and dynamic equilibria are 
compared for complex fractions of 50 and 80 percent. The dif- 
ferences are so huge that even without fitting the curves we 
can note the main important outcomes: 1) Considering a static 
equilibrium, the SACFs (gray and light magenta curves) have 
only diffusion terms whereas the SACFs for a dynamic equilibri- 
um (green and blue) decay much faster due to an additional 
bunching term; 2) the average number of molecules per spe- 
cies is equal in both cases (static and dynamic equilibrium) but 
all molecules in the mixture are contributing to the diffusion 
term of the dynamic equilibrium. This is nicely demonstrated 
for the case of 50% percent complex (Figure 4 B) by overlaying 
all SACFs, where the data of the complex (light magenta 
squares) hide those of the ligand (gray squares). Each of them 
contains the same particle number in the diffusion term as de- 
fined in Equation (18) [1/(A/(1-X (i) )) = 20]. We fitted the SACFs 
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with a model function similar in structure to one usually used 
for fitting the diffusion with a single bunching term [Eq. (18)]: 



steady state concentrations of species which we can determine 
in the next step of our analysis (see Section 2.2.2). 
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where N is the total number of labeled particles in the excita- 
tion volume, X® is the fraction of species /, and f R is the relaxa- 
tion time of binding equilibrium (the other parameters are as 
defined before). 

In this model there are no contributions defined by changes 
in brightness like the singlet-triplet transition case where the 
molecule is either bright or dark. In this case, the molecule cor- 
responds to species 1 or 2. The fits show the feasibility of this 
function. The obtained amplitudes agree nicely (error < 3%) 
with simulations, where the total number of molecules in the 
focus contributing to diffusion was 0.1 for both species. In the 
case of 80% complex (Figure 4C) the amplitudes of the corre- 
lation curves for both species are now clearly different and 
closely reflect the values expected from the simulations. 

To check the existence of a ligand-receptor dynamic binding 
equilibrium, we used SCCF. The SCCFs for three bound frac- 
tions (20%, 50% and 80% complex) are presented in Fig- 
ure 4D. The existence of an exchange dynamics is nicely high- 
lighted by the presence of an amplitude of SCCF (colored 
lines), whereas the amplitude is zero for a static mixture (gray 
line). The SCCFs were fitted using the following model function 
assuming equal diffusion times for all species (f d = f d '' = f d m) ) 
[Eq.(19)]: 



G^(t c ) = 1 +lj-G d „(tc) ■ (1 - exp(-f e /f»)). 



(19) 



In Equation (19), N is the total number of labeled particles in 
the excitation volume and f R is the relaxation time of the bind- 
ing-unbinding process (the other parameters are defined as 
before). The overlay of SCCF and fit curves immediately shows 
that SCCF alone cannot distinguish between 20/80% and 80/ 
20% mixtures, since SCCFs have equal amplitudes and relaxa- 
tion times. But, of course, these two cases have different 



(18) 



.kjj 2.2.2. Global Analysis of SACF and SCCF 

The SACF contains interrelated information on 
species concentrations and relaxation time 
whereas the SCCF contains additional selective 
information on the relaxation time. The kinetic 
term of SACF is scaled by the equilibrium con- 
stant, where it is not in the SCCF; that is, if the 
equilibrium is shifted to one side, the amplitude 
of the kinetic term of the two SACFs will vary 
characteristically, whereas it won't change in the SCCF. There- 
fore global analysis of SACFs and SCCF will stabilize the fits sig- 
nificantly. If we express our model functions in terms of rate 
constants directly, global target fit of SCCF and SACFs will in- 
crease the accuracy of analysis and immediately will provide 
the values of rate constants. This is of particular interest for 
studying binding processes. To accomplish global target analy- 
sis we have to combine Equations (18) with (19) and replace f R 
by (^1,2 + ^2,1) Finally, we get the following model functions 
[Eqs. (20)]: 

= 1 +i-G diff (0- (1 -exp(-(/c,, 2 +/c 2i1 ) -f c )) 



^2,1 



exp(-(/c 1-2 + /c 2J ) -t c ) 



d L '\t c ) = 1+^-G dm (tc)-(l + 
G m M = 1 + lj ■ G m (t c ) ■ (1 + ^ ■ exp(-(/c,, 2 + k v ) ■ t c )\ 



(20) 



The results of global target fits for the example of the simulat- 
ed dynamic binding equilibrium of a ligand as described in 
Section 2.2.1 (Figues 3 and 4) are compiled in Table 1. The rela- 
tive deviation between simulated and fitted values ranges be- 
tween 0.1 and 3%. Given the fact that the obtained relative 
deviations are determined by the variations of simulated 
values and the confidence intervals of the fits, we computed 
by sampling of data subsets 1501 and extrapolation to larger 
photon numbers a theoretically expected error of less than 
2%. The excellent agreement with the simulated parameters in 
Table 1 confirms the high accuracy and precision of the ap- 
plied analysis. 



Table 1. Global target fit results of SACF and SCCFs for 50% 


and 80% formed complex 


in simulated dynamic binding equilibria (Figures 2 and 3). 


Parameter 




50% bound 






80% bound 






Simulated 


Fit results 


Rel. deviation [%] 


Simulated 


Fit results 


Rel. deviation [%] 






Free ligand Complex 






Free ligand Complex 




,V FCS 


0.05; 0.05 


0.0513 0.0507 


2.6; 1.4 


0.02; 0.08 


0.0204 0.0797 


2.0; 0.4 


f d [ms] 


1.0 


1.021 


2.1 


1.0 


1.002 


0.2 


k u [ms- 1 ] 


2.0 


2.011 


0.6 


8.0 


8.057 


0.7 


k v [ms- 1 ] 


2.0 


2.034 


1.7 


2.0 


2.059 


3.0 


W FCS 


0.1 


0.102 


2.0 


0.1 


0.1001 


0.1 
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2.2.3. Temporal Resolution of Polarization-Resolved fFCS to 
Detect the Conformational Exchange between Two States 

To explore the temporal resolution of polarization-resolved 
fFCS, we simulated typical polarization-resolved fluorescence 
decays of a single fluorescent dye attached to a macromolecule 
(Figure 5 A). The time-resolved anisotropy r(f) is described by 
three exponents [Eq. (21)]: 



r(t) = Vv exp(-f//o overa n) + r 0 .ba exp(-f/p backbone ) 
+r 0 ,u exp(-t/p 

linker/ 



(21) 



with the characteristic times for overall rotation (p overa n), back- 
bone fluctuations of the macromolecule (p backbone ) ar, d rota- 
tions around the linker (p| inker ), respectively, fractional initial ani- 
sotropy amplitudes r 0J . The steady-state anisotropy r ss is com- 
puted by Equation (22): 
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Figure 5. A) Conditional probabilities of polarization-resolved fluorescence 
decays used for fFCS. Total brightness Q tot = 400 kHz, r = 4 ns, W = 0.1, para- 
meters of the observation volume are as in Section 2.2.1. The details on r(f) 
are given in the text. B) The species cross correlation functions [data (black 
circles: f R = 100 us, open black circles: f R = 10 us, black squares: f R = 1 us, 
open black squares: f R = 0.1 us), fits by Equation (19) (blue lines)] nicely re- 
cover the simulated parameters (magenta). 



We simulated a case where the anisotropy senses the confor- 
mational dynamics of macromolecule as observed for a flexible 
loop in bacteriorhodopsin by Schroder etal. [511 In the simula- 
tion, we studied a change of p backbone due to the assumed con- 
formational exchange in macromolecule, that is, state 7 charac- 
terized by a flexible backbone is in equilibrium with state 2 
characterized by rigid backbone. In state 2 the effective frac- 
tion of the overall rotation is actually represented by the sum 
of fractional initial anisotropy amplitudes (r 0iOV + r 0iba ). In the 
simulation we assumed a fluorescence lifetime r = 4 ns and the 
following rotational correlation times p, (with corresponding 
fractional anisotropy amplitudes r 0J ): p overaN = 20 ns (0.18); 
Pbackbone 01 =1 ns (0.10) and Pbackbo„e (2) = 20ns (0.10) for state 
1 and 2, respectively; pi inker = 0.3 ns (0.10). The conformational 
change results in a typical increase of the steady state aniso- 
tropy from r ss (1) = 0.177 to rj 2) = 0.240. To find the temporal res- 
olution limit we reduced the relaxation time for the exchange 
by factors of ten starting from 100 us and ending at 100 ns. If 
the properties of the two states are known, the filter functions 
for conditional probabilities can be directly used to calculate 
the corresponding SCCFs shown in Figure 5 B. Analysing ex- 
change between the two states by Equation (19), the simulat- 
ed dynamics can be easily recovered over at least three orders 
of magnitude reaching a sub-microsecond time resolution. 

Similar anisotropy differences in r 5S are also observed if flexi- 
bly linked dyes stick temporally to the macromolecule, which 
may lead to complications in a quantitative analysis of FRET ex- 
periments. 1521 fFCS and anisotropy Photon Distribution Analysis 
(PDA) 1531 open up the possibility to detect heterogeneous dye 
environments and to characterize their exchange in an experi- 
ment, so that these results can be compared with detailed MD 
simulations. 



2.3. Experimental Results for Single-Molecule FRET 
Experiments 

One exciting experimental application of fFCS in biophysics is 
the study of protein conformational dynamics. As an example, 
we have analyzed Syntaxin 1 (Sx); Sx forms part of the so- 
called SNARE (soluble NSF attachment receptor) proteins that 
regulate synaptic vesicle release. 154,551 Sx is a transmembrane 
protein and the soluble domain consists of four long alpha 
helices that are known to undergo a conformational transition 
upon interacting with targets. One of the targets is SNAP25 lo- 
cated in presynaptic vesicles and it is known to open Syntaxin 
1 allowing vesicle fusion and neurotransmitter release there- 
after. On the other hand, Munc-18 is expected to act as a nega- 
tive regulator of exocytosis. 155,561 However, Munc-18 is also es- 
sential for exocytosis and maybe it also works as an activator 
instead of inhibitor. 1571 Even when free in solution Sx had previ- 
ously shown to be in equilibrium between two conforma- 
tions, 1341 in the following referred to as "open" and "closed". Sx 
has two domains, the Habc domain consists of a three-helix 
bundle and a flexible H3 domain has a single long helix. H3 as- 
sociates with the Habc domain in the closed and dissociates in 
the open conformation. 
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Herein, we used fFCS to ex- 
tract the kinetic rates of this 
transition in a single-molecule 
FRET experiment combining life- 
time and polarization-resolved 
information to maximize the 
contrast between conforma- 
tions. The Sx double mutant 
(G105C/S225C) was labeled with 
donor dye Alexa488 and accept- 
or dye Alexa594 as described in 
Section 4.4.' 341 A solution of 
double-labeled protein was di- 
luted to -10 pM in PBS and 
placed on the MFD setup (Sec- 
tion 4.3). After careful calibration 
of the detection efficiencies of 
our experimental setup, we per- 
formed an analysis of the single- 
molecule bursts at first. The re- 
sults are displayed in two-di- 
mensional frequency histograms 
of the FRET indicator F D /F k (ratio 
of donor fluorescence over ac- 
ceptor fluorescence) 1491 against 
fluorescence lifetime of the 
donor in the presence of accept- 
or, r D(A) . In Figure 6A, one can 
see a smeared population cover- 
ing a broad range of fluores- 
cence lifetimes. At least two 
FRET populations in addition to 
the third small D-Only contribu- 
tion at high F D /F A ratios with t d(0) 
of 4 ns can be detected (Fig- 
ure 6 A). Such a presentation 
allows one to check for devia- 
tions from the static FRET line 
(orange line). The static FRET 
line in F D /F A versus r D(A) plots is 
defined as [Eq. (23)]: 
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Figure 6. A) 2D Histogram F D /F A versus r D(A1 of Sx 105/225 labeled with Alexa488 and Alexa594. Grayscale contours 
represent binned accumulation of single-molecule events with the following parameters: The average green and 
red background count rates were (6 g ) =4.4 kHz, (S,) =0.39 kHz, respectively. An estimated 4.2% of cross-talk 
signal was accounted for. The green over red detection efficiency (g G /9n) was n -8 and the green quantum yield 
<?> FD(0) = 0.8, red quantum yield 0 FA = O.4, D-Only lifetime r D(0) = 4.0 ns, and g-factors 0.995 and 1.378 for green and 
red channels respectively. On the histogram the main FRET population clearly shows a distribution that is off the 
static FRET line given by Equation (23) (orange). In dashed green, the dynamic FRET line Equation (24) shows the 
path taken by the conformational exchange. In this case there are two conformational states, plus a small fraction 
of D-Only. The open state was identified to have a lifetime T D(A) lop<in ' = 3.6 ns and rotational correlation time of 
p= 1.5 ns. For the closed state r D(A1 lclo! ' 6dl = 0.8 ns. B) Conditional probabilities used for generating the filters: buffer 
(IRF), open [r D(A) lop,inl = 3.6 ns,p = 1.5 ns] and closed [r Dlfl) (dosedl = 0.8 ns,p=1.5 ns], C) Species cross-correlation be- 
tween closed and open states. Data is fit with Equation (19) and shows a relaxation time of f R = 0.6 ms, and 
f d = 1.8 ms. On top of this panel, residuals of fit with one or two relaxation times are shown. Additionally, we fit 
this SCCF globally with SACFs, shown in D), and the fit required two relaxation terms instead of one (f R , = 1.1 ms, 
f R2 = 0.08 ms) with the same f d = 1.8 ms. D) SACFs of the open and closed states as defined by the filters described 
above. Fits represent the global target fit of the two SACFs and the SCCF. The difference between global target fit 
of the SACF and SCCF with one and two relaxation times are shown on top of the SACF. The fit requires two relax- 
ation times to reduce the characteristic deviations in the residuals. The SACFs fit with one relaxation time are not 
shown. 



where r D(A) determined by the maximum likelihood estima- 
tor/ 591 corresponds to the fluorescence-weighted average life- 
time {t d(a) = (r) f ). For further details on the analysis and inter- 
pretation of MFD histograms, we refer the reader to Sisamakis 
et al. [491 The smeared population clearly follows the dynamic 
FRET line shown in green. The end points of the dynamic line 



(23) 



where <P Fm and <t> l 



represent the quantum yields of the 
donor and acceptor fluorophores, respectively. The fluores- 
cence lifetime of the donor without acceptor is t d(0) and r DW 
with an acceptor. As the populations are clearly off the static 
line, dynamics faster than the diffusion time is present. 1581 Then, 
given the simplest case of a dynamic two-state system, a FRET 
line for this exchange (green dashed line) can be traced that 
follows the form [Eq. (24)]: 
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correspond to r D(A) (open) = 3.6 ns and r D(A) (dosed) = 0.8 ns, which 
were determined by sub-ensemble analysis of the dynamic 
FRET population. In combination with the average rotational 
correlation (p=1.5 ns, also reported in ref. [34]) we computed 
two fluorescence decays with Equation (37) and the corre- 
sponding conditional probabilities (Figure 6 B), which were 
used to generate the filters according to Equation (11) using 
the global error minimization methodology described in Sec- 
tion 2.1. The SCCF computed by Equation (12) is presented in 
Figure 6C. One can clearly see the enhanced negative correla- 
tion term as in the case of simulation data. When fitted with 
Equation (19) one recovers a relaxation time f R = 0.6 ms which 
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agrees nicely with values of f R = 0.8 ms [341 and f R = 0.6 ms [491 re- 
ported previously for another Sx double mutant (S91C/S225C), 
which also probes the exchange dynamics of the H3 domain 
with respect to the Habc domain. 

In the previous section, we advised to use global target fit 
of the SACF and SCCF to stabilize the fit and reduce uncertain- 
ties. As indicated by the residuals of SACFs in Figure 6 D, the 
global target fit of Sx revealed the need of an additional relax- 
ation time found in all three correlation functions simultane- 
ously with following relaxation times and formal amplitudes: 
f Rl = 1.1 ms (84%), f R2 = 0.08 ms (16%). The f R = 0.6 ms obtained 
with one free parameter recovers the average behavior, but 
only the global fit shows the need of the second faster relaxa- 
tion time. The f R2 = 0.08 ms term may be caused by a more 
complex kinetic scheme for the conformational transition. For 
example, Reiner et al. [601 have postulated that large-scale con- 
formational motions require an intermediate structural unlock- 
ing step leading to a more flexible state which reacts to the 
final open state. Further measurements are needed to unam- 
biguously solve this question. Finally, it is worth mentioning 
that the temporal resolution of fFCS is so powerful and pro- 
vides such a good contrast that it is possible to distinguish ad- 
ditional kinetics not possible to observe otherwise, even with 
the 2D histograms from smFRET data. 

If we assume an effective two state model for the open- 
closed transition, we can combine the results of the SACF and 
SCCF and extract the kinetic rates. Using the target fit Equa- 
tion (20) with the average relaxation time of 0.6 ms, the indi- 
vidual rate constants /c, 2 = 0.62 ms" 1 and k v =0.94 ms -1 and 
the corresponding equilibrium constant K 0/c ='\.5 (W (open) = 0.09, 

/V (dosed) = 0 06) a| . e obtainec | 

To conclude, Sx was selected as an experimental realization 
where FRET adds the additional spectral and lifetime informa- 
tion to maximize contrast in fFCS. The combination of the 
SACFs and SCCF has shown to be a powerful tool to extract 
the kinetic rates and the equilibrium constant between con- 
formers. 

3. Discussion 

The introduction of lifetime filtering opened a door to increase 
contrast on FCS curves by selectively interrogating species 
rather than photons. In this paper we extend this idea to in- 
clude polarization information to increase sensitivity. One 
could add also spectral information for selectivity as it will be 
discussed. During these studies several differences to previous 
implementations' 22 " 24,611 were observed and are worth discus- 
sing. 

3.1. Species Distinction by MFD: Generation of the Most 
Efficient Filters 

One example is the simulation done in Section 2.1. In this case, 
the fluorescence of two species with the same fluorescence 
lifetime but small difference in their rotational correlation 
times (p (,) = 0.1 ns and p <2) = 0.3 ns) is studied. It was now the 
task of finding the most efficient filter. Scenario 1: single detec- 
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tor, single filter {' d ~ 5 f): If the signal is detected by a single detec- 
tor (see Figures 1 A and 9 A, B) and the excitation is polarized, 
the polarization induced difference in the fluorescence decays 
is very small, so that conventional fluorescence lifetime correla- 
tion spectroscopy hardly accomplishes to separate the species 
(in Figure 1 A). Scenario 2: two detectors, independent filters 
( 2d "'"r): The separation is improved by using polarized excitation 
and polarization-resolved detection; one detector for the paral- 
lel and one for the perpendicular component (see Figures 1 B 
and 9C, D). If the error minimization for the filters is done for 
each channel independently the separation of species im- 
proves, but the information on the relative intensities of polar- 
ized signals is lost. Scenario 3: two detectors, global filter 
( 2d " 9 'f): Using polarized excitation and polarization-resolved de- 
tection, most efficient filters are generated by global error min- 
imization (see Figures 1 C and 9E, F) by properly maintaining 
the relationships between the parallel and perpendicular de- 
tection channels. 

One way of judging the efficiency of the filters is to compare 
the quality of the computed SACFs by quantifying the noise 
by average normalized / 2 values with respect to the model 
functions. If we set the % 2 values of the most efficient method 
with global filter ( 2d_9 'f) equal to one, the noise increases by 
a factor of 6.3 for independent filters ( 2d_ "Y) and even by 
a factor of 200 for a single filter ( ld ~ s f). Moreover, we checked 
the capability of the filters to separate species by calculating 
the particle number N and the species brightness Q = F/N. The 
best agreement between simulated and obtained values was 
found for the global filter ( 2d ~ 9 'f). To conclude, all three ap- 
proaches for fFCS work in principle but filters, which use the 
resolved signal from multiple detectors and are generated in 
a global way, are by far the most efficient. 

In the following we discuss what features of the underlying 
species-specific fluorescence decays make the resulting filters 
most efficient for fFCS. The relative difference of the fluores- 
cence parameters (e.g. r or p) forms the basis to distinguish 
two species. Considering the currently used organic fluoro- 
phores and detectors in confocal microscopes, the fluores- 
cence lifetime can range only between 0.2 and 6 ns, whereas 
the dynamic range of the rotational correlation time is signifi- 
cantly larger (0.2 and at least 30 ns), which offers a large dy- 
namic range for the contrast. Furthermore, if we consider dif- 
ferences solely in the rotational correlation time, this difference 
should ideally affect the whole decay curve and not only 
a short segment as, for example, shown in Figure 7A. In other 
words, if p m and p (2) r, the two decays are nearly equal for 
any time t>p (Figure 9E), no matter how different the two ro- 
tational correlation times are. This results in equal conditional 
probabilities to detect a photon characteristic for a certain spe- 
cies and thus limits the discrimination efficiency. We demon- 
strate this effect by comparing two simulations, where the ro- 
tational correlation time of species 1 is kept constant (p (1) = 
0.1 ns) and the rotational correlation time of species 2 is either 
p (2) = 0.3ns (Figure 9E) or p (2) =1.0ns (Figure 7A). The low 
noise of SACF in Figure 7B with p (2) =1.0 ns is consistent with 
drop of the average normalized Rvalue by a factor of 100 as 
compared to SACF with p (2) = 0.3 ns (Figure 1 C). 
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Figure 7. fFCS separates species when difference in rotational correlation 
time between two species is 0.9 ns whereas the fluorescence lifetimes are 
the same (4 ns). All other simulation conditions are the same as in Figure 1. 
The total duration of the simulated experiment is 1598 s. A) Conditional 
probability of the two species [species 1 p <1) = 0.1 ns in blue and species 2 
p (2> = 1.0 ns in red], IRF shape is shown in gray for reference. B) fFCS can sep- 
arate species. SACF of species 1 and 2 are shown in blue and red respective- 
ly. The simulated species are shown as dashed lines for the same case, and 
represent the expected correlation functions. The noise is significantly re- 
duced from the case presented in Figure 1. 



In cases where the differences of a chosen parameter (e.g. 
rotational correlation or lifetime) are not large enough, the ad- 
dition of another experimental observable can be used to in- 
crease contrast. For example, one can combine changes in life- 
times, rotational correlation times and spectral information 
(like FRET) to be more precise in the selection. 1331 In this case 
up to four decay curves (two for donor and acceptor, respec- 
tively) can be combined in a global filter to maximize contrast 
and thereby minimize noise. 

3.2. Influence of the Filter Quality 

To make a precise selection, there are some experimental re- 
quirements and limitations that are worth discussing. It is well 
known 1451 that the response function of avalanche photodiodes 
(APDs) exhibits a noticeable dependence on the photon count 
rate, which can distort the decay patterns for individual spe- 



cies. This could lead to small time shifts on the measured fluo- 
rescence decay of the mixture compared to the decay of indi- 
vidual species. We observe this shift in measurements with 
a mean particle number N = 0A, because under these condi- 
tions the count rate is much higher than under the conditions 
where we measure the IRF by using the Raman and Raleigh 
scattered photons of the pure buffer solution. In Figure 8 we 
show two SCCF and their corresponding fit curves to demon- 
strate how the shift by 0.3 ns and 0.2 ns in parallel and perpen- 
dicular detection channel, respectively, can affect the baseline 
in the correlation curve. At shorter correlation times the SCCF 
curve with shifted decay patterns enters the range of negative 
correlation amplitudes which makes no scientific sense. How- 
ever, the relaxation times in these two cases are not affected. 
This outcome is quite positive, because it means that if we 
allow the fit to enter into negative correlation amplitudes 
range we can still recover the relaxation rate between the two 
species. 




Correlation Time / ms 

Figure 8. Overlay of SCCFs calculated by fFCS for dynamic mixture of species 
with X"- R) = 0.5 [free ligand with p lftee " gan<l) =2.71 ns and complex with 
^(complex) _ 5 1 g ns j in bu ff er w jth out (light green) and with shift (dark green) 
of species decay curves relative to IRF. Further simulation conditions are 
given in Figure 3. Shifts of patterns are 0.3 ns and 0.2 ns for the detection 
channels of parallel and perpendicular polarization, respectively. 



Another problem is to find appropriate filters that describe 
each separated species. In many cases the properties of the in- 
dividual states can be determined by shifting the equilibrium 
towards one side, for example by adding a binding partner or 
stabilizing/destabilizing the macromolecule. In the case of slow 
dynamics it might be possible to resolve the species in 2D 
plots and generate sub-ensemble decay histograms. 1611 More- 
over, a set of ensemble measurements with high photon statis- 
tics will also help to judge the heterogeneity of the sample. 
These results can be then compared with other methods such 
as 2D plots and PDA which are very useful to detect structural 
dynamics by FRET. 1581 If slightly wrong decay curves are used to 
generate filters for fFCS, the amplitudes of SACF and SCCF 
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become distorted but so far we can still recover correct relaxa- 
tion times. 



n{=2) 



(2) 



1=1 

n(=2) 



(25) 



3.3 Strategy to Eliminate Detector Afterpulsing Effects in 
Mixtures 

Another experimental artifact is the occurrence of highly corre- 
lated events referred to as afterpulsing of the detectors up to 
the time range of 10 us. It has been shown by Enderlein 1231 that 
FLCS is able to correct the correlation curves of a single spe- 
cies for afterpulsing in a single detection channel by adding an 
uniform baseline as a model for "signal" generated by after- 
pulsing. In our experience, this technique does not work well 
for the case of multiple species mixture. Thus, SCCF and SACF 
from real data will be affected by afterpulsing in the short cor- 
relation times range (less than 10 us). According to ref. [62] we 
can avoid the afterpulsing correlation and detector dead time 
effects by splitting the fluorescence signal between two inde- 
pendent detectors and cross-correlating these two signals. This 
requires the addition of a 50/50 beam splitter in our parallel 
and perpendicular detection channels. The theory described in 
Section 2.1 is the same, but we need to correlate the signals of 
one parallel/perpendicular detector pair with the correspond- 
ing signal from second parallel/perpendicular detector pair 
using their corresponding filter functions for each pair of de- 
tectors. In total four green detectors, instead of two green de- 
tectors in a typical MFD setup, have to be used. With this solu- 
tion one could increase the temporal resolution to below us. 
This implementation and the exploitation of this powerful 
technique are left for further studies. 

As a final remark we want to give a recommendation for op- 
timal measurement conditions of fFCS. Due to the filter action, 
the effective ("useful") signal is significantly lower than the raw 
signal so that the noise levels of all filtered FCS curves for 
given experimental conditions are significantly higher than 
that of conventional FCS curves. If the average number of mol- 
ecules in observation volume is chosen too high, the noise 
level of the SACFs and SCCF becomes comparable to the am- 
plitudes of diffusion and kinetic terms of the correlation func- 
tion and the analysis becomes ambiguous. This problem can 
be avoided by using concentrations, which are optimal for 
fFCS measurements with an average number of molecules N< 
0.1, which was also chosen in simulations. This concentration is 
slightly lower than for standard FCS, because the background 
signal is anyhow filtered out. 



where "p y w and ± pj' ) are the normalized probabilities of count- 
ing a photon within the y'th TCSPC channel for species / in par- 
allel (| |) and perpendicular (±) polarizations relative to linearly 
polarized excitation laser beam respectively. To compute the 
normalized probabilities we follow ref. [22] for a single detec- 
tion channel where the counts in the TCSPC histograms {H^} 
were normalized by the total number of photons [Eq. (26)]: 



Pj 



E H f 



£/f = i 



i-1 



(26) 



with ;'=1, 2 ; y = 1 ;. . ; L; and L is the number of histogram bins 
or TCSPC channels. In the case of two detectors, one parallel 
and one perpendicular, the independently normalized proba- 



bilities, 



and ± pf> 



E ll p (0 



are [Eq. (27)]: 



Pi 



(27) 



In Equation (25) w l,) is the number of photons correspond- 
ing to the rth species. To increase selectivity, Equation (25) 
could be expanded to account for spectral differences (color 
detection) by including a summation over the number of color 
channels. 

Two filter sets and L fj, with /=1; 2, can be derived to 
satisfy the relationships [Eqs. (28)]: 



H; 



Hw(') 



£i- w 



H; 



i-1 



(28) 



and x w w are the photon counts of the fth species 
in parallel and perpendicular detection channels, respectively 



w w ). Brackets represent averaging over an in- 



finite ensemble of measurements. The relative errors in both 
detection channels are simultaneously independently mini- 
mized like [Eq. (29)]: 



4. Theory and Methods 

4.1 Theory for Independent Error Minimization (Scenario 2 
of Filter Generation) 

TCSPC decay histograms {"H y } and { x Hj} of a mixture of two 
species (/=1; 2) can be written as a superposition of species 
histograms in the form of [Eq. (25)]: 



L 



min and 



(29) 



then ^fj' and ± ff' > are expressed with the help of the weighted 
pseudoinverse of the matrices "M and 



Ma 



L pf ), as 1431 [Eq. (30)]: 



M, ("M„ = "pf and 
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H/f> = (["M-cliag^H,)- 1 - "M 7 | "/Of • diag("tfy)" 



x f = [ X /W ■ dlag( x H y )- 1 - x Jtf r ] - X M ■ diag( x H ; 



(30) 



since single-photon detection in a TCSPC channel follows Pois- 
sonian statistics. In Equation (30), diag("Hy) _1 and diag( x H,) _1 
are L x /.-dimensional diagonal matrices with diagonal ele- 
ments ("H,-)" 1 and ( x Hy)~\ 7=1;...; L, and 7 represents the 
transpose operator. It can be shown that x fy W and the pat- 
terns "pj , x p|' are orthonormal systems [Eq. (31)]: 



1 , same species : i = k, 

0, distinct species : /^/c. 

1, / = *, 

0, 



(31) 



With the filters ^f™ and x ry (,) it is now possible to separate the 
correlation curves of two molecular species independently. The 
fFCS functions between perpendicular and parallel fluores- 
cence signals G]*Jt c ) are defined by Equation (32): 



E x f . x Sy(f)J ■ [py-hj(t + t c ) 



(32) 



where x S,(f) is the perpendicular signal in the y'th TCSPC chan- 
nel at measurement time f and "5y(f + f c ) is the parallel signal 
in the yth TCSPC channel at measurement time f + f c . The po- 
larization cross-correlation functions of a single species 
/ G^||(f c ) and G[| ± (f c ), averaged over an infinite number of 
measurements or over sufficiently long time {>t c ), can be writ- 
ten as [Eq. (33)]: 



r 1 " 



( x F(--)( f) -^)(r + fc )) 
{cl ( x Fm(f))-(iiF«(f + t c )) 



("F")(t) ■ x F"-»(f + f c )) 
{iiF(0(t))-( x F(0(r + O) 



(33) 



where x F (,, (r) and "F (,) (f + r c ) are pure fluorescence signals 
from the /th-type molecular species (summed over all TCSPC 
channels). The global fit of two correlation curves significantly 
reduces the statistical errors and makes fit results more accu- 
rate. 

To show how to calculate the filters according to the theory 
presented we use the simulated mixture introduced already in 
Section 2.1. A mixture of two fluorescent species with the 
same (4 ns) lifetime, but different characteristic rotational cor- 
relation times (0.1 ns and 0.3 ns) in water are simulated. 

Given polarized excitation and single detection channel at 
magic angle conditions (54.7°) separation of species with the 
same lifetime is impossible. Even in MFD, using the total decay 
histograms as h t = "f/y + 2 X F/, (equivalent to magic angle con- 



ditions) does not provide separation of species. However, if un- 
polarized detection (assuming h t =" = x /?y = "Hy + x H y ) is 
used, one would expect small differences between TCSPC 
decay histograms of two species. Figure 9A, B shows the step- 
wise generation of the filters for this case imitating single de- 
tection channel experiment. 

First, the conditional probabilities are required for each spe- 
cies [Eq. (26)]. These can be generated independently by meas- 
uring fluorescence decays (Figures 9A,C,F) by TCSPC with pure 
concentrated solutions of each species (>25 nM; on average 
about 30 molecules in a 2 fl detection volume). However, if the 
species are not available separately, it is also possible to com- 
pute the theoretical decay patterns based on the expected 
fluorescence lifetimes and rotational correlation times. This can 
be accomplished by measuring the instrument response func- 
tion (IRF), convoluting it with the expected fluorescence decay 
for each species and normalizing these histograms using Equa- 
tions (4)-(6), (26), (27). The single detection channel filters are 
then obtained by minimizing the errors according to Equa- 
tions (11), (29). The blue line corresponds to species 1 with 
p n) = 0.1 ns and the red line corresponds to species 2 with 
p (2) = 0.3 ns. Finally, the SACFs [Eq. (12)] for the case of a single 
detection channel, with no split in polarization, using the filters 
shown in Figure 9 B, are calculated (Figure 1 A). There is poor 
separation of species and their correlations (blue and red lines 
for the species p (,) = 0.1 ns and p (2) = 0.3 ns respectively) are 
very noisy. 

The separation of species gets better when the two polar- 
ized decay components are considered as shown in Figur- 
es 9C,E. In Figure 9C, we can see already the advantage of 
measuring the distinct orientations. The normalized TCSPC 
decays are distinct due to 0.2 ns difference in the rotational 
correlation times between species. For two independent detec- 
tion channels the error minimization described by Equa- 
tion (29) is performed. These filters are then used to calculate 
the SACFs for species by Equation (12). The independent filters 
also were stacked and applied to the raw data in a modified 
format where the TCSPC channels of both detectors are 
stacked in a single array with length 2L. This correlation curves 
show better separation of species than in the single detector 
case. 

Figure 9E shows conditional probabilities for each species as 
described in Equations (4)-(6). The IRF shape is shown in gray. 
The first species has a rotational correlation of p (l) = 0.1 ns and 
its decay is shown in blue. The second species with rotational 
correlation p <2) = 0.3 ns is shown in red. The filters generated 
for each species (fj w ) according to Equations (10), (11) are plot- 
ted in Figure 9 F. The IRF filter is shown in gray, species 1 in 
blue, and species 2 in red. Calculated by Equation (12) SACFs 
are shown in solid blue and red lines for species 1 and 2, re- 
spectively. For comparison, the simulated correlation functions 
are shown as dashed lines in dark blue and wine for species 1 
and 2. The separation of species is obvious to the eye, and the 
estimated error in the recovered parameters is within 2%. In 
all cases, we make sure that exactly the same number of pho- 
tons are used and the quality of separation is defined only by 
the different filters. 
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Figure 9. Decays and filters for simulated mixture of two species with same fluorescence lifetime but different ro- 
tational correlation times. The parameters used for each species in the simulation are Q (,I = Q I2) = 150 kHz, 
t d m = t d m = 0.25 ms, N l1l = 0.133, N l2l = 0.266, dark counts = 0.2 kHz, scatter = 2.5 kHz. The lifetime of both species is 
the same (t g "' = t g (21 = 4.0 ns), and the respective rotational correlation times are p (,1 = 0.1 ns, p (2 ' = 0.3 ns. The total 
duration of the simulated experiment is 2479 s. A) The normalized TCSPC histograms for') for each species as de- 
scribed in Equation (26). The IRF shape is shown in gray. The first species (blue line) has p' 11 = 0.1 ns and the 
second species (red line) has p 12 ' = 0.3 ns. B) The generated filters for each species (f] 1 ') are plotted according to 
the notation: IRF shown in gray, species 1 in blue, and species 2 in red. C)The independently normalized TCSPC 
histograms for species 1 and 2 according to Equation (27). D) Simultaneously independently minimized filters 
using Equations (29) and (30). E) Conditional probabilities for each species as described in Equations (4)-(6). The 
IRF shape is shown in gray. The first species has a rotational correlation of p" ' = 0.1 ns and its decay is shown in 
blue. The second species with rotational correlation p l2l = 0.3 ns is shown in red. F) The filters generated for each 
species (fP) according to Equations (10) and (11) are plotted. The IRF filter is shown in gray, species 1 in blue, and 
species 2 in red. 



proved for a given box size; ii) a 
fractional number of molecules 
can be realized, which is impor- 
tant for irreversible reactions; 
and iii) the simulation speed can 
be improved because the calcu- 
lations can be easily parallelized. 

In our simulation of an "open" 
volume, molecules are allowed 
to leave the simulation volume, 
after which they are not further 
tracked. To keep the desired 
average number of molecules 
constant, new molecules are 
added at each time step as de- 
scribed below. 

Analog to the 3D Gaussian ob- 
servation volume in Equation (3), 
the simulation volume is defined 
by a surface l(x,y,z) = const, 
which for our observation 
volume can be approximated by 
an ellipsoid [Eq. (34)]: 



m 



(34) 



with m = 10 unless stated other- 
wise. One can show that the 
average number of molecules 
entering the simulation volume 
per time step Af is given by 
Equation (35): 



(N ln ) = cSy/D d At/jt 



(35) 



while the distribution of distan- 
ces A/ of new molecules from 
the ellipsoid surface Equa- 
tion (34) is [Eq. (36)]: 



Al 



d/ 



(36) 



4.2 Data Generator 

Our simulations were done by the Brownian dynamics ap- 
proach 163 " 661 with an own in-house C+ + program. The spatial 
intensity distribution of the observation volume was assumed 
to be a 3D Gaussian. 1671 In previous implementations, Brownian 
diffusion was modeled in a box with periodic boundary condi- 
tions. 1631 We have developed our own approach, where "mole- 
cules" diffuse freely in an "open" volume. In this way, we gain 
several advantages: i) the efficiency of the algorithm is im- 



In Equations (35) and (36), 
c stands for the concentration of the molecules, S is the sur- 
face area of the ellipsoid Equation (34), and D d stands for the 
diffusion coefficient. Equations (35), (36) can be obtained by 
solving the diffusion equation with initial conditions c(/<0, 
t=0) = c and c(/>0, f=0) = 0. The curvature of the surface Equa- 
tion (34) is thereby neglected, which requires that the time 
step Af is small (i.e. D d Af ma> 0 ). The number of new mole- 
cules per time step (Af = 0.005 ms) was generated according 
to a Poisson distribution with a mean value given by Equa- 
tion (35). Test simulations have shown that the number of mol- 
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ecules calculated by FCS analysis Equation (1) of simulated 
data agree with the assumed values within 0.2-0.5%. Note 
also that the average number of molecules in the simulation 
volume can be fractional, which is not the case for a box with 
periodic boundary conditions. As a result, all concentration 
values are possible irrespective of the simulation volume size. 

Our approach allows also for straightforward and unbiased 
modeling of photobleaching and other irreversible reactions: 
when bleached molecules leave the simulation volume, new 
"bright" molecules are automatically added, Equations (35) and 
(36). Finally, the trajectories of molecules in "open volume" are 
not endless, which allows efficient distribution of calculations 
between several CPUs. 

To model diffusion, normally distributed random numbers 
Ax, Ay and Az were added to the x, y and z coordinates of 
each molecule, respectively, at each time step Af, where 
(Ax 2 ) = (Ay 2 ) = (Az 2 ) = 2D d At. The average number of photons 
emitted by a molecule and registered by the /th detector at 
each time step is (F,) = Q, • /(x,y,z) ■ Af, where Q, stands for 
the brightness of the molecule as measured in the /th detec- 
tion channel. 

Binding or conformational dynamics were modeled by al- 
lowing states 1^2 and 2^1 transitions. The initial fractions of 
molecules in states 1 and 2 were ^2,1/(^1,2+^2,1) anc ' K.il 
(/c 12 +/c 21 ), respectively. The times the molecules spent in 
states 1 and 2 (f, and f 2 , respectively) were exponentially dis- 
tributed with P(f 1 ) = /c 12 _1 exp(-/c 12 f 1 ) and P(f 2 ) = 
/c 21 _1 exp(— /c 21 f 2 ). The diffusion times of the species are given 
for each case and the time step was 0.005 ms for all simula- 
tions. Background was added when necessary. 

TCSPC data were generated as proposed by Chowdhury et 
al., [681 assuming a Gaussian excitation pulse profile with an 
FWHM of 0.3 ns. The parallel and perpendicular decay compo- 
nents F| |(f) and F ± (f) were modeled as [Eq. (37)]: 

F || (t)=F(f)-(1+(2-3-/,)T(t))/3 
F ± (f)=F(f).(1-(1-3./ 2 ).f(f))/3 

In Equation (37), F(f) is the fluorescence intensity decay (typi- 
cally mono-exponential with r = 4ns), /, = 0.0308 and / 2 = 
0.0368 are correction factors, 148,691 and the anisotropy decay r(f) 
is either single-exponential or given by Equation (21). The 
background signal consists of dark counts (uniformly distribut- 
ed over the TCSPC channels) and a scatter contribution [re- 
peating instrument response function (IRF) of the setup]. The 
simulated data were saved in the data format of SPC-132 
TCSPC cards (Becker & Hickl GmbH, Berlin, Germany) and were 
suitable for direct software testing. The algorithm was thor- 
oughly tested by analyzing the simulated data by the FCS, 
FIDA, and b I FL [9, 16,45,641 methods. 

4.3. MFD Experiments 

The MFD setup used for single-molecule FRET detection was 
previously described. 1321 In short, the active mode-locked 
argon-ion laser was operated at 477 nm (28 kWcrrr 2 ). A dual- 
band dichroic mirror 488/594PC (AHF Analysentechnik, Tubin- 



gen, Germany) was used to separate excitation light from fluo- 
rescence. The detection spectral windows for the donor (D) 
and the acceptor (A) were defined by bandpass filters HQ535/ 
50 and HQ650/75 (AHF Analysentechnik) respectively. Burst in- 
tegrated lifetime analysis using the signal of two green and 
two red detectors split by polarization were detected by ava- 
lanche photodiodes (AQ 151, EG&G (Vaudreuil, Quebec, 
Canada)) coupled to a counting board (SPC432, Becker & Hickl, 
Berlin) and a personal computer. Single-molecules bursts were 
selected, analyzed and used to generate the corresponding 
sub-ensemble TCSPC histograms for the mixture and the MFD 
2D histograms. 

4.4. Molecules 

Syntaxin-1 was expressed from pET28a in the BL21 (DE3) strain 
of Escherichia coli as previously reported. 170,711 The recombinant 
Sx contained N-terminal His6 tags that served to affinity purify 
it on nickel-nitrilotriacetic acid-agarose (Qiagen). After affinity 
purification, DTT (5 itim) and EDTA (2 him) were added. The 
tags were cleaved off by thrombin. Cleavage occurred over- 
night during concomitant dialysis (NaCI (50mM), Tris HCI 
(20 mM, pH 7.4) and DTT (1 itim)). Subsequently, Sx was puri- 
fied with MonoQ column using the Akta Explorer (Amersham 
Pharmacia Biotech). Collected fractions were analyzed by SDS/ 
PAGE and Coomassie staining with 95% purity. Before labeling, 
DTT was removed from the proteins by gel filtration using PD- 
10 columns and PBS as elution buffer (Na 2 HP0 4 (10.4 itim), 
KH 2 P0 4 (3.2 mM) and NaCI (123 mM)). 

Random labeling of Sx was accomplished by adding a mix- 
ture of Alexa 488- and Alexa 594-maleimide dyes (Molecular 
Probes) with a ratio of 10:30 vs. protein concentration respec- 
tively in order to obtain an 1:1 labeling. After 3-5 hr of reac- 
tion on ice, labeling was stopped by adding DTT (10 mM). The 
unreacted dye was removed via gel filtration using a PD-10, 
followed by extensive dialysis against PBS (Na 2 HP0 4 (10.4 mM), 
KH 2 P0 4 (3.2 mM) and NaCI (123 mM)) containing DTT (1 mM) 
and SM-2 Biobeads (Bio-Rad). All measurements were done in 
PBS buffer previously cleaned by charcoal. 

5. Conclusions 

fFCS extends the previously presented FLCS 122 " 24,611 by including 
the anisotropy and spectral information for improved species 
separation. With some modification, fFCS is capable of discrimi- 
nating detector afterpulsing from fluorescence signal as has 
been described in Section 3 as well as removing effect of the 
scatter photons on the correlation amplitude (Section 2). More 
importantly, we have shown that fFCS is able to distinguish 
species even when they have very close or even the same fluo- 
rescence lifetime using just differences in time-resolved fluo- 
rescence anisotropy, with an estimated error of <2%. In addi- 
tion, to increase selectivity, one can combine parameters such 
as lifetime, spectral window and anisotropy. One example is 
the protein Synatxin 1, which shows conformational dynamics 
in solution between open and closed conformation. A relaxa- 
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tion time of 0.6 ms recovered by fFCS agrees with previous re- 
sults. 

One has to bear in mind that there are different error mini- 
mizations forms and this is totally different from the previously 
presented lifetime FCS. [221 Given the particular case, it is impor- 
tant to consider the differences on minimizations procedures. 

Finally, standard tools such as FCS and FCCS could not even 
distinguish the static binding from dynamic one if there is no 
brightness change upon binding, as shown in Section 2, and 
fFCS now opens a window to study kinetic mechanisms in 
detail and we foresee many possible applications: 1) Removal 
of background: In experiments at very low molecular concen- 
trations and for measurements in highly scattering media, fFCS 
can be used as tool to eliminate the influence of a scattered 
background on the FCS curves. Usually, the drop of amplitude 
is corrected by using a formula of Koppel. 1141 However, we have 
shown this correction works only to a limited extent/ 161 where- 
as fFCS fully recovers the amplitudes, which is crucial for 
a proper concentration determination. Also in FRET-FCS experi- 
ments with a dynamic equilibrium the correct extraction of 
amplitudes or species concentrations is the prerequisite to de- 
termine the equilibrium constant. 2) Selective correlation in mix- 
ture: fFCS can be applied as a tool to compute species-specific 
correlation curves from mixture of different species, provided 
that a separation could be done based on lifetime or/and ani- 
sotropy or/and spectrum. 3) Analysis of interconversion kinetics 
between species: In our group, SCCF and SACF have developed 
as a hallmark for kinetics, which also allows for the assignment 
of the originating states. The applications can be split into two 
general groups: i) Experiments with single-color detection 
based on anisotropy changes due to ligand binding (static or 
dynamic), homo-FRET and mobility changes of a dye coupled 
to a biomolecule (e.g. dye sticking observed in FRET studies), 
ii) Experiments with multicolor excitation and detection based 
on lifetime and anisotropy changes of species (in all single- 
molecule FRET studies of biomolecules with MFD technique). 
In this case, differences in the FRET efficiencies of the species 
resulting in different green and red signal levels make it possi- 
ble to include the spectral heterogeneity as a third parameter 
for species selection. 
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